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Abstract. Spectroscopic and imaging data of the peculiar galaxy Mkn 298 are presented in this paper. Narrow-band H a and 
broad-band R images are used to study the star formation rate in the galaxy and its morphology, which is typical of a merging 
system. Long-slit and integral field spectra are used to assess the kinematics of gas and stars, and the nature of the ionizing 
source at different distances from the nucleus. In particular, the nucleus of Mkn 298 is characterized by peculiar line ratios: 
[Nn]/i6583/H a is typical of Hn-like regions, while [OfJ/16300/Ha' could indicate the presence of an active galactic nucleus. 
We show that models where a shock component is added to photoionization from a starburst allow to reproduce the observed 
line ratios. Mkn 298 is thus most likely a star forming galaxy, rather than a galaxy hosting an active nucleus. 
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1. Introduction 

Mkn 298 (IC 1182) is known as a morphologically pecu- 
liar system (Figure located in the Hercules supercluster. It 
shows a chai n (c, d, e) of small compact blue emitting re- 
gions dStocktonl 19721) aligned on its eastern side up to ~ 80" 
from the main body of the galaxy, which is classified in NED 
(NASA/IPAC Extragalactic Database) as a SA0+ galaxy. These 
regions are aligned with a tail (b) resembling a spiral arm lo- 
cated on the eastern side of the galaxy (a). On the western side 
a very faint trace of spiral arm is also visible. A f aint emis- 
sion is visible in the c knot, b ehind a foreground star ( Stockton 
Il972l lMetik &Proniklll982t) . The d and e knots are charac- 
terized by an emiss ion line spectrum. They were identified by 
Brai ne et all feOOlb as a tidal dwarf galaxy (TDG hereafter) 
forming from material ejected from the disk of a galaxy after a 
collision. 
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* Based on data obtained at the 2.2m and NTT telescopes of ESO-la 
Silla (Chile), at the 2.2m telescope of DSAZ-Calar Alto (Spain), and 
at the 6m telescope of SAO (Russia). 



The classification of Mkn 298 as a gala xy hos ti ng an 
acti ve nucleus is controversial in literature. iKoskil dl978l) 
and IVeilleux & Osterbrockl Jl987l) classified it as a Seyfert 
2 galaxy with an extremely low [Nnl/16583/Ha line ratio. 
Conversely Viegas-Ald rovandi & Gruenwa ld ( 1990) enclosed 
Mkn 298 in a sample of LINERs. Variability was claimed 
bv lMetik & Pronikl Jl982t) who observed changes of the rel- 
ative intensity of [Om]/Ll4959,5007 and HyS in the period 
1969 -70, 1977. The d etection of variable soft X-ray emission 
(Rafan elli etai"1ll999h with a luminosity typical of Seyfert 2 
galaxies, seems to co nfirm the presence of an active nucleus. 
However, as noted bv lMoles et ai1 ( l2004 it is not possible to 
exclude that t he X-ray em ission is produced in star-fo rming re- 
gions dZezas et al. 120021) . Finally. Moles et alJ [2004) detected 
the presence of several knots in the nucleus; they interpreted 
Mkn 298 as a very luminous starburst (SFR ~ 90 MQ/yr) with 
very low metallicity (Z/Zq ~ 0.06 - 0.1). 

An extended Hi distribution pea ked on th e TDG 
was detected usin g VLA data jBraine et all l200lt 
llglesias-Paramo et alJ 120031) . iBraine et alJ J200ll) report 
the detection of CO(2-l) and CO(l-O) emission in Mkn 298, 
however the emission is peaked on the main galaxy and is 
below detection in the TDG. They computed a total amount of 
H2 mass in the system M mo \ > 5 x 10 9 Mq, and an upper limit 
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M m0 ] < 6 x 10 7 Mq in the TDG. This would imply that only 
a small fraction of HI was transformed into H2 in the TDG, 
while in the nucleus the conversion already took place. This is 
consistent with Mkn 298 being in the late stage of a merging: 
as outlined by iHibbard & van Gorkoml ([1996), in advanced 
mergers the atomic gas is mostly relegated to the outer tidal 
features and is almost absent in the central regions, whereas 
the molecular gas is concentrated in the remnant body. 

The aim of this paper is to study in detail the physical con- 
ditions in the different components (galaxy body and tail) of the 
Mkn 298 system by the analysis of their emission-line spectra. 
Observations and data reduction steps are described in Sec. [3] 
The morphology, star formation rates and kinematics are ana- 
lyzed and discussed in Sec.|3j Sec. |4] and Sec. [5] respectively. 
In Section [6] emission line ratios are used to derive chemical 
abundances. Observed ratios are then compared with the re- 
sults from both pure photoionization models and photoioniza- 
tion+shock models, using the Cloudy and Suma codes respec- 
tively. Conclusions are finally drawn in Sec. [71 

2. Observations and data reduction 

Emission line spectra and images of Mkn 298 have been taken 
in different periods and using different observing facilities. 

2.1. Imaging 

1. Narrow band images in He (/I e =6820 A, A/l=170 A) and 
continuum light (i c =6240 A, AA=l 10 A) were taken at the 
Deutsch-Spanische Astronomische Zentrum of Calar Alto 
(DSAZ) 2.2 m telescope (Spain) in September 1996, with 
a 2kx2k CCD (spatial scale 0.53"/px), and under a seeing 
of ~ 1.5". Images of the galaxy of 600 and 900 s were 
obtained with the two filters, respectively. Moreover, the 
spectrophotometric standard star LDS749b was observed 
for flux calibration. All images were processed by means 
of the IRAF 1 packages. They were firstly bias subtracted, 
flat-fielded and cosmic -rays cleaned. Then, we performed 
the photometry of the standard star, whose fluxes were mea- 
sured with PHOT, corrected for atmospheric extinction, and 
compared to the fluxes obtained from the spectrum of the 
star convolved with the response curve of the filters, in or- 
der to obtain the calibration constants. The filter transmis- 
sions were conveniently taken into account. The images of 
the galaxy were aligned, background subtracted, normal- 
ized to 1 s exposure and corrected for atmospheric extinc- 
tion. Then, the continuum image was scaled by a factor 1 .6 
to compensate the different shapes of the two filters. The 
respective calibration constants were applied to the so pro- 
cessed images, and finally the continuum was subtracted 
from H a to obtain a pure emission line image. 

2. Five broad-band R consecutive exposures of Mkn 298, 15 
min each, were obtained with EMMI at the ESO-NTT 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



in August 1991. The 2kx2k CCD, with 15yU-size pixels, 
yielded a spatial scale of 0.348"/px. The seeing was 1.4". 
The five images were bias and flat-field corrected, then 
aligned and combined to obtain an average frame cleaned 
from cosmic rays. No photometric calibration was per- 
formed. 

3. A public image of the galaxy was extracted from the HST 
Archive. Obtained in June 1995 with the WFPC2 in combi- 
nation with the F606W filter, and 500 sec of exposure time 
(see iMalkan et al.lll998i) . this image was simply cleaned 
from cosmic-rays and directly analyzed. The main part of 
the galaxy was covered with the PCI, having a spatial scale 
of 0.045 "/px. 

2.2. Spectroscopy 

1. The nucleus was observed in May and June 1999 with 
the Multi Pupil Fiber Spectrograph (MPFS) at the 6 m 
telescope of the Special Astrophysical Observatory (SAO 
RAS, Russia). The MPFS is an integral field spectrograph, 
in which an array (16x15 at the epoch of the observations) 
of l"xl"microlenses is connected to a bundle of fiber op- 
tics. The 240 fibers are re-arranged to form a pseudo-slit 
which is the input to the spectrograph. A 600/mm grat- 
ing was used during the first run (Rl) of observations, in 
combination with a lkxlk CCD. Each of the two spectral 
ranges, 3350-6100 A and 4550-7300 A, was observed with 
a dispersion of ~ 2.6 A/px, an instrumental FWHM of 7.5 
A (~ 400 kms~'at 5500 A) and an exposure time of 3600 
sec. During the second run (R2) a 1200/mm grating was 
used, which covered the ranges 4800-6200 A and 5900- 
7300 A, with a dispersion of ~ 1.36 A/px and an instrumen- 
tal FWHM of 4.5 A (~ 250 km s 1 at 5500 A). Integration 
times of 7200 sec for the "blue spectrum" and 4800 sec for 
the "red spectrum" were applied. The seeing was typically 
around 1.5-2" for both runs. 

The MPFS spectra were reduced and flux calibrated with a 
special software developed in IDL at SAO. 

2. Knots d,e were also observed with the MPFS during the 
above runs. In this case, a 300/mm grating was used, which 
allowed to cover the whole optical spectral range 3600- 
9000 A, with an instrumental FWHM ~ 15 A. An expo- 
sure time of 3600 sec was applied. Long-slit spectra of the 
same knots were taken with EFOSC2 at the 2.2-m telescope 
at ESO-La SiUa (Chile) in April 1995. The instrumental 
FWHM is ~ 10 A, the covered range is 4500 - 7000 A. 
The standard data reduction steps (bias subtraction, flat- 
field correction and flux calibration) were applied using 
the LONGSLIT package in IRAF. Due to the higher spec- 
tral resolution of the ESO spectra, we decided to use them 
rather than the MPFS data in the analysis of the knots, with 
the exception of the [On]/l3727 line which is covered by 
the MPFS spectra only. 

The positions of the slits and MPFS spectra are displayed 
in Figs.^ an d|2] where they are overplotted to the R-band and 
H a images of Mkn 298 respectively. The H a image clearly 
shows the existence of two knots in the nucleus which were 
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Fig. 1. The Mkn 298 complex in the /?-band image with over- 
laid the position of the MPFS arrays. 

fu rther resolved in fo ur knots (CK,CKN1 and CKN2,CKN4) 
bv lMoles et all fe004 . The knots are not visible in the 7?-band 
image. 

A 2" x 2" rebinning of the MPFS spectra was done before 
the measurements in order to increase the signal to noise ratio: 
each measured spectrum is therefore the sum of four MPFS 
spectra. Spectra from R2 were used to obtain a 2D map of 
the emission line ratios with the highest spectral resolution. 
Finally, two regions corresponding to the nuclear knots (Nl, 
N2) were defined and spectra inside them were summed to- 
gether. In this case, spectra from Rl were used to cover the full 
spectral range from [On]/l3727 to [S ii]/16716, 31, but with a 
lower spectral resolution. 

Emission line ratios after corrections for reddening are 
shown in Tabled For each region, we selected the spectra with 
the best spectral resolution, namely the MPFS spectra for the 
nucleus and the ESO spectra for the d,e knots. Uncertainties on 
the measured line ratios are < 10% in the nucleus, ~ 20-30% in 
the knots. The uncertainty is much higher (~ 50%) in the case 
of [Om]i4363, due both to the nearby Hy /14340 and to the 
noise in the continuum. 

The redshift of Mkn 298 as derived from the peak positions 
of the lines is z — 0.034, giving a linear scale 0.66 kpc/arcsec 
(H = 75 km s" 1 Mpc" 1 ). 



3. Morphology 

Understanding the peculiar morphology of Mkn 298 requires a 
detailed analysis of deep images. For this purpo se we app lied 
to the NTT R-band image the GIM2D package (iSimard et alJ 
2002j), which performs a two-dimensional bulge/disk decom- 
position of galaxy images, and produces as output the model 
convolved with the PSF and the residual image (= galaxy - 
model). Mkn 298 was first fitted using a Sersic /?'^"-law for the 
bulge and a classical exponential law for the disk. The fit gave 
as result an index n = 3.5 and a reduced-^ 2 ~ 1.28. Then a de- 
Vaucouleurs 7?'^ 4 -law was applied for the bulge obtaining quite 
similar residuals (reduced-^ 2 ~ 1.16). The resulting B/T ratio 
~ 0.78 indicates clearly that Mkn 298 is a spheroidal galaxy, 
likely a SO, with a weak disk inclined at i ~ 40°. The fit gave 
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Fig. 2. Continuum-subtracted Ha image of Mkn 298 with 
overlaid the position of the ESO slit on the d,e knots. 



# 




Fig. 3. Residuals after the GIM2D two-dimensional model sub- 
traction: white arrows indicate the long tail (LT), the dust lane 
(DL), the plume (P) and the shells (Sh). 



also the effective radius of the bulge, r e ~ 1" and the radius 
Hmif ~ 9.5" containing half of the total integrated light. 

The structure and morpholog y of the ionized gas in the nu- 
cleus was discussed in detail by Mo les et all J2004h . In addi- 
tion, we note that the strong nuclear dust lane, whose maxi- 
mum radial extension is ~ 17" (~ 11 kpc), leaves a residual 
resembling a ring inclined at i ~ 140 - 145°. Interestingly, 
the galaxy shows at least three shells located at ~ 40", 1' and 
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Fig. 4. HST image of Mkn 298 and fit of the surface brightness 
profile. 



1.3'(~ 26, 42 and 50 kpc respectively) at East, North-East to 
the nucleus, and another one at ~ 58" (~ 38 kpc) North- West to 
the nucleus, connected to a large plume o riented at P. A. ~ 160°, 
and named Second Tail bv lMolesetdl fc004 (see Figure |3j- 
These features are not uncommon in elliptical and SO galax- 
ies, and are still studied to understand their real natu re and ori- 
gin. The simulations bv lHernauist & Soergell l l 19921) show that 
mergers between equal-mass gala xies can deve l op su ch shells 
on timescales of about 1-2 Gyr. Turnbull et al. ( 1999) studied 
two examples of shell ellipticals, which have very strong analo- 
gies with those present in Mkn 298 (see their Figs. 9 and 10). 
Nevertheless, Mkn 298 is a rare example of early-type galaxy 
showing both shells and a prominent stellar and gaseous tidal 
tail, a situation not seen in the above me ntioned numeric al sim- 
ulations. Faint tidal tails were found by Balcells ( 1997) in the 
shell elliptical NGC 3656, and several shells and a long H i 
tidal tail were discovered bv IVerdes-M onteneg ro et all y002) 
in HCG 54. Only the NW shell has a surface brightness close 
to 3cr above the background. Assuming a sky brightness ~ 
20.5 mag arcsec -2 , that is the value expected at LaSilla for 
similar moon illumination conditions, we obtain fj. R ~ 25 mag 
arcsec~ 2 : such a value is in agreement with those measured 
by Turnbull et al. for the shells of their galaxies. The other 
shells in Mkn 298 have surface brightness between 1 and 2<x 
of the background. In the residual image it can be also seen 
that the Eastern tidal tail extends beyond the d and e knots, up 
to 1.75'(~ 70 kpc), where it seems to turn back toward North- 
West. 

The nucleus of Mkn 298 was also observed by HS T with 
the W FPC2 and the broad-band V filter F606W jMalkan et alJ 
Il998l) . An analysis of the isophotes of the galaxy was done 
using the ELLIPSE task in IRAF. Since the presence of the 
dust lane makes difficult to fit the real isophotes, we applied 



Table 1. Ha luminosities and star formati on ra t es. Th e labels 
in parentheses are those adopted by Mo les et alJ |2004). 



Observed Dereddened 





Log L Ha 


ceo 


Log L Ha 


ceo 


Nl (CK+CKN1) 


41.46 


2.30 


41.83 


5.35 


N2 (CKN2+CKN4) 


41.21 


1.29 


41.50 


2.52 


a (galaxy body) 


41.72 


4.19 


42.01 


8.14 


b (MTK1) 


40.02 


0.08 






c (MTK3) 


< 39.87 


<0.06 






dl (MTK4) 


39.93 


0.07 


40.52 


0.26 


d2 (MTK5) 


40.03 


0.08 


40.58 


0.30 


e (MTK6) 


40.42 


0.21 


40.42 


0.21 



the masking method used by ICarollo et al The sur- 

face brightness profile plotted against r 1 ^ 4 can be fitted by a 
de Vaucouleurs law (r e ~ 8") without an additional disk com- 
ponent; it does not show any nuclear point-like source (e.g. 
AGN or stellar cluster), in agreement with the results given by 
GIM2D. 

In summary, the brightest part of Mkn 298, that is within 
~ 6.3 kpc, is clearly spheroidal. Such morphology, combined 
with the presence of a dust lane, a bright and extended tail and 
some weak shells, are all strong clues in favor of the hypothesis 
that Mkn 298 has experienced a major merger event. 

4. Star formation rates 

The presence of multiple knots in the nucleus is also visible 
in the continuum subtracted H a image. Aperture photometry 
was done on the two brightest nuclear knots (Nl, N2) and on 
the regions in the tail (b,c,d,e). A very faint Ha emission is 
detected in the c knot. 

An estimate of the photometric accuracy was done by com- 
paring the so-obtained fluxes with those obtained from the 
MPFS spectra. We conclude that the photometry of imaging 
and spectroscopic data is consistent within 5%. 

Ha luminosities obtained in each of these regions are 
shown in Tabled for those regions where dust extinction from 
Ha/HB was derived from spectra, the extinction corrected 
values are also displayed. Ha luminosities were f inally con- 
verted to star formation rates using Equation 2 from lKennicutt I 
(1998). In the nucleus, this is an upper limit if part of the H a 
emission is due to an AGN. We see that the total star formation 
rate (a region) is ~ 8.1 Mp yr" 1 , we ll below the value of 90 
Mq yr~' given by Mol es et al.l(l2004J). The Ha lumin osities in 
the galaxy measured by us and bv Moles et al 1 ll2004 before 
correction for dust reddening are consistent within 30%. The 
difference i n the SFR is mainl y due to the internal reddening 
adopted by Moles et all d2004 . E(B - V) ~ 1, which is much 
higher than the average reddening that we find in the nucleus, 
E(B — V) ~ 0.3. This value of the r eddening is also in agree- 
ment with the value found bv lKoskildl978l) . E(B - V) = 0.40. 
A possib le reason for this discrepancy may be that the spec- 
tra that iMoles et"afl d2004l) used to derive the extinction come 
from knots (CK and CKN2) which are very close to the dust 
lane: they could be therefore biased to high values of H a/HB, 
not representative of the whole galaxy. 
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Fig. 5. Velocity fields of Ha (up-left), [Om]/i5007 (up-right), and stars (bottom-left); the cross indicates the position of the 
nuclear continuum. North is up and East on the left, (bottom-right) The deprojected velocity curves of Ho- (filled circles), 
[Om]/i5007 (open squares) and stars (open stars). The solid line is the pure circular motion fit of the gaseous component. 



Mkn 298 is not a bright IRAS source; the IRAS Faint 
Source Catalog contains an entry which is however centered 
on the d,e knots. Approximate fluxes at the position of Mkn 
298 can be derived from the Iras Sky Survey Atlas using the 
IRS A 2 facility: we obtain f^,,, = 0.8 Jy, /ioo//m =2.1 Jy. These 
values are very uncertain as the emission is distributed on a 
wide area (> 2 arcmin around the galaxy) and is probably the 
sum of the emission from the galaxy and from the tidal tail. 
According to Equations 7, 8 from Hopki ns et al.l (120031) . we 
obtain L(8 - 1000/zm) = 1.8 x 10 44 erg s' 1 and a SFR ~ 1 1 M© 
yr _1 , a value in good agreement with the SFR derived from the 
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total Ha luminosity (galaxy +tidal tail). As the contribute to 
the SFR from the tidal tail estimated from the H a luminosity 
is negligible compared to that in the galaxy, it seems likely that 
this is also the case for the infrared luminosity. 

iGeorgakakis et alJ lkoOO) explored the possibility to use the 
L F i R /M(H2) ratio (star formation efficiency, SFE hereafter) as 
an indicator of the galaxy merging age. In the case of Mkn 
298, we obtain SFE < 9 L©/M© using the above far-infr ared 
lumin osity and the value M mo \ ~ 5 x 10 9 M© from Br aine et"all 
d200 ll) . Together with the low fraction of HI in the nucleus, this 
value would indicate that Mkn 298 is in a post-merger stage. 
The Nl and N2 knots in the nucleus are probably separate star 
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forming regions partially obscured by the dust lane, rather than 
separate nuclei before merging. 

5. Kinematics 

Taking advantage of the integral field data, we carried out a 
two-dimensional investigation of the gaseous and stellar kine- 
matics in Mkn 298. 

First, the wavelength positions of the two brightest emis- 
sion lines [Om]/15007 and Ha were measured by means of 
gaussian fitting of their profiles in each spectrum where they 
could be detected down to 3cr over the continuum. Then, the 
Mgb A5175 absorption triplet was chosen for deriving stel- 
lar kinematics, and the Fourier cross-correlation metho d was 
applie d by means of the IRAF task FXCOR ( Ton rv & Davisl 
1 19791 TD79 hereafter), and using a mean sky-fiat spectrum as 
template. Only the correlation peaks with the TD79 parameter 
R > 5, which is an indicator of the signal-to-noise ratio, were 
considered reliable. The so obtained radial velocities were cor- 
rected for the motion of the observer in the direction of the ob- 
servation, and then re-organized to construct the velocity fields 
for each component (Figs. [5}- A first look reveals that the [O m] 
and H a kinematics are in agreement, showing similar displace- 
ment of their "spider" diagrams, which are extended up to a 
maximum radius of ~ 7" (~ 4.6 kpc), and concordant veloc- 
ity values. In addition, their velocity fields have the major axes 
clearly parallel to the dust lane, strongly suggesting that the gas 



and the dust lane belong to the same structure, a rotating disk 
of a likely former spiral galaxy merged into a early-type, or 
however a poor-gas galaxy. On the contrary, stars show a ve- 
locity field with a slightly different orientation, indicating the 
presence of a kinematic al decoupling. 

In order to determine the rotation curve of the galaxy in 
the field of view of the spectrograph, we applied the GIPSY 
(Groningen Image Processin g SYstem) task R OTCUR. Based 
on the method described by BegemarJ dl989h . ROTCUR de- 
rives kinematical parameters by fitting the velocity field with 
concentric tilted-rings, characterized by seven free parameters. 
In our case, given the relatively small number of available mea- 
surements, it was necessary to convenientl y reduce the deg rees 
of freedom. We followed the recipes by Begeman ( 1983) to 
determine the Xq and Yq coordinates of the center, and the 
systemic velocity, V jys ~ 10250 kms~', and in addition we 
fixed the inclination, i, of the rings on the basis of the optical 
isophotes. Therefore, ROTCUR fitted the circular velocity, V c , 
for increasing radii, varying the position angles of the rings. 
We performed this procedure separately for the approaching 
and the receding sides of the velocity field, in order to ac- 
count for possible distortions or asymmetries. The resulting 
V C (R) values were plotted to visualize the deprojected rota- 
tion curve for gas and stars. As mentioned above, [Oin] and 
H a show similar radial velocities, almost symmetrically dis- 
tributed, with a AVmax ~ 250 kms -1 . Stars are well in agree- 
ment with the gaseous component only in the receding side, 
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Table 2. Dereddened emission line ratios (H j3- 1 .00) and models computed assuming (a) photoionization only and (b) photoion- 
ization + shock (see text for model parameters). The values in italics are the uncertainties on the observed line ratios; the H a/H f3 
ratio before reddening correction is also given below. 





MPFS 
Nl 


MPFS 
N2 


ESO 
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ESO 
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Cloudy 
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AGN + shock 
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SB + shock 
Nl N2 d,e 
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1 from MPFS data 



while poorly in the approaching side, with lower velocities up 
to -100 km s~' with respect to the gas. From the analysis of the 
position angles of the rings, we obtained that the kinematical 
axes of gas and stars are not aligned, as expected, and oriented 
at ~ 30° and ~ 0° respectively. Moreover, the kinematical axis 
of the gas is more stable in the receding side than in the ap- 
proaching one, where the variations can reach 20°. 

We have com pared our rotation curve to that published by 
iMoles et"al] l|2Q04), finding a general good agreement for the 
gaseous component. These authors used a long-slit spectrum 
taken close to the major axis, and obtained a curve with similar 
radial extension and velocities. Thanks to their higher spatial 
resolution they could claim about two rotating systems, which 
unfortunately are not visible with our data. On the contrary, we 
can assert that the stellar component shows rotation, at least in 
the region of the galaxy where the stellar kinematics could be 
measured. 

We fitted the H a curve with an analytic relation, which as- 
sumes a spherical potential model with pure circular motions, 
V C (R) = aR/(R 2 + c 2 Q y /2 jBertola et alll99ll) . 

The nonlinear least-squares fit was stable, giving the fol- 
lowing parameters: a = 515 kms -1 , co = 3.8"and p = 1.27. 
Using this relation, we estimated the total mass within a radius 
of 7"obtaining ~ 7.3 x 10 10 M©. 



The nuclear stellar velocity dispersion, cr , was also mea- 
sured. First, the brightest central spectrum was considered, and 
then an average of central spectra within an aperture of 3". 
The results were similar. The estimate of <tq was obtained 
again with FXCOR by measuring the FWHM of the correla- 
tio n peak, and then convert ing into <x following the recipes 
bv lNelson & Whittle! J 1995b . In detail, we convolved the tem- 
plate spectrum with gaussian functions of increasing known 
<T, and we cross-correlated each convolved spectrum with the 
original one, obtaining a sort of cr-FWHM empirical relation, 
which becomes flat for <x < 90 km s _1 . The nuclear cr resulted 
~ 170 kms" 1 , a va lue in agreement with the one measured by 
Mol es et alJ (l2004). even if slightly lower. 

Given the existence of a Mg-<x relation in early-type 
galaxies and spiral bulges, we measured the LICK/IDS 
system Mg and Fe indices in the center o f the galaxy: 
Mgi,M g ? , Mg b, Fe521Q, Fe5335, defined by IWorthev et all 
(1994), and we attempted a comparison with published val- 
ues and relations. Since emission lines are sometimes present 
at the blue or red edges of these absorption features, the es- 
timate of the continuum may be rather uncertain. Therefore, 
our measurements should be considered as rough estimates. 
Anyway, Mg\ = 0.08, Mg 2 = 0.17, Mg b = 3.6 and 
< Fe >= 0.5(Fe5270 + Fe5335) = 2.07 are values typi- 
cally observed in spiral bulges, and not in ellipticals or SOs 
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Fig. 7. VO diagnostic line ratios in Mkn 298. Circles: binned 
MPFS spectra in the nucleus. Dots: integrated MPFS spectra in 
the nucleus (Nl, N2). Triangles: ESO spectra in the d,e knots. 
The solid line is the theoretical starburst border from Kewley 
et al. (2001), dashed lines are +0.1 dex of these lines. 

(see e.g. lProctor & Sansomll2002l) at the ctq measured in Mkn 
298. In other words, Mkn 298 is an outlier of the cr-Mg and 
cr-Fe relations, showing indices with lower strength with re- 
spect to spheroidal galaxies. This effect can be caused either 
by age or by metallicity. On the other hand, the [Mg/Fe] ra- 
tio is around 0.4, indicating that the bulk of the nuclear stellar 
component is old probably formed in a timescale of 1-2 Gyr 
ftdiartet.al]ll99fil> . We can conclude that Mkn 298 does not 
fit the properties of the early-type galaxies because of the on- 
going merger. Indeed, bursts of star formation induced by the 
strong interaction can effectively make the Mg and Fe indices 
strengths weaker because o f the presence of bright, hot stars 
JWorthev & Colloberfcool) . 

Similar results were obtained by Longh etti et alJ |2000), 
who showed that shell- galaxies, believed to be the product of 
a past strong interaction, deviate significantly from the cr-Mg 
and cr-Fe relations (see their Fig. 2). They proposed that this 
may be due either to the gravitational interaction which affects 
the velocity dispersion or to bursts of star formation that alter 
the metallicity indices, or to a combination of both. 

6. Emission line ratios 

Figure0shows the VO line ratios observed both in the nucleus 
and in the d,e knots: overlaid are the lines dividing the theo- 
retic al_sterburs^egiOTifrom other types of excitation computed 
bv iKewlev et al.l (12001 ). While in the d,e knots the line ratios 
are consistent with those expected in a starburst, the situation 
in the nucleus is more complex. In fact, in both the nuclear 
knots [Nn] /I6583/Ha < 0.3 is typical of starburst regions, 
[Sn] AA 6716,6731 ~ 0.5 may be fitted by starburst models, 
whereas [Oil/Ha > 0.1 is too high for any starburst model. 
A similar behavior may be seen in all the individual spectra 
taken by the MPFS array in the nucleus. Figure[8]and Table|3] 



Table 3. Statistics of dereddened emission line ratios in the 
nucleus from the MPFS data (Rl). 





min 


max 


avg 


a 


[On] ,13727 


3.53 


9.78 


7.39 


1.89 


[Om] ,15007 


0.78 


2.90 


1.69 


0.61 


[Oi] ,16300 


0.27 


0.63 


0.40 


0.09 


Ho- ,16563 


2.18 


2.80 


2.76 


0.13 


[Nn] ,16583 


0.47 


0.99 


0.72 


0.14 


[Sn] ,16716 


0.56 


1.08 


0.79 


0.12 


[Sn] ,16731 


0.41 


0.63 


0.51 


0.06 


[Sn] ,167 16+31 


1.01 


1.65 


1.30 


0.16 


E(B-V) 


0.03 


1.58 


0.50 


0.42 



show that line ratios are remarkably constant over the whole 
galaxy (up to 5" from the center). Reddening by dust is low 
(E(B — V) < 0.50), with the exception of the regions on the 
western side (E{B - V) > 1), where it is due to the foreground 
dust lane. 

6.1. Chemical abundances 

As discussed in more detail in the next section, emission line 
ratios in the nucleus may imply the presence of an active nu- 
cleus. Chemical abundances in the nuclear regions were there- 
fore computed using both the two calibrations proposed by 
IStorchi-Bergmann et alJ Jl998i SSCK her eafter) for t h e NLR 
of active galaxies and the R23 metho d dMcGaughl Jl99ll) 
iKobulnickv et al.l J 1 9991) ^ and p method JPilvuginl2000l l2001) 
for abundances in H 11 regions of normal galaxies. We did not 
use [Om]/14363 due to the strong uncertainty in its measure- 
ment. 

The SSCK calibrations are based on [Nn]/U6548,6583/Ha 
vs. [Om]AA 4959,5007/H/? (Eq. 2 in SSCK) and 
[On]A3727/[Om]AA 4959,5007 (Eq. 3 in SSCK). Following 
the SSCK prescription we adopted the average of the 
two values, 12 + log(0/H) ~ 8.7 (Z/Zq ~ 0.7, assuming 
12+log(0/H)Q = 8.87). To check t hese results we use d the 
intensities given in Table 5 from iMoles et"al] J2004I) and 
obtained consistent results, 12 + log(0/H) ~ 8.5 (Z/Zq ~ 0.6). 

The /?23 and p methods allow to compute oxygen abun- 
dances from [Om]AA 4959,5007/[On]/l3727. Both give as re- 
sult a lower abundance, 12 + log (Q/H) ~ 8.3 (Z/Z n ~ .3). 

As it concerns the d,e knots, ll glesias-Paramo et alJ d2003h 
found 12 + log(0/H) = 8.41 (Z/Z ~ 0.35) using the p-method. 
A similar value, 12 + log(0/H) = 8.38 (Z/Zq ~ 0.33) is ob- 
tained from the emission line ratios that we measure in the e 
knot. 

We therefore conclude that there is evidence for a sub-solar 
metallicity throughout the galaxy. The metallicity in the nu- 
cleus (Z/Zq ~ 0.5-0.7) may be somewhat higher than in the d,e 
knots (Z/Zq ~ 0.3). We find no evidence supporting t he ver y 
low metallicity (Z/Zq ~ 0.06) claimed by M oles et 2112004). 

Starting from these chemical abundances, we explored dif- 
ferent models trying to fit the observed line spectra: pure pho- 
toionization models relative to an AGN and to a starburst, and 
composite models which account for both photoionization and 
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shocks. In all the following models, Z/Zq was initially set to 
the value found here and then varied to give the best fit of the 
observed line ratios. The model which leads to the best fit of 
the line ratios will determine the controversial nature of Mkn 
298. 

6.2. Photoionization models 

As a first attempt to model the emission line ratios measured 
in the different regions, we assumed that they are produced 
by photoionization from either an active nucleus or stars (H n 
regions). To this ai m, the photoionization code Cloudy 94 
jFerland et al J 19981) was used. 

In order to find the best-fit parameters from photoioniza- 
tion models, we run a grid of models with different values of 
ionization parameter and hydrogen density. 

An attempt to use as the ionizing source a Kurucz contin- 
uum with temperatures up to 50000 K confirmed that in this 
case it is not possible to fit the [O i]/H a ratio. We then took as 
the ionizing source a power-law continuum (table power law 
option in Cloudy), F v oc y~ L5 . 

Starting from the solar abundances given in CLOUDY 3 , 
abundances of all elements heavier than helium were scaled 
as Z/Zq, with the exception of nitrogen which was scaled as 
(Z/Z Q ) 2 . 

We considered acceptable a model if line ratios (H/3-l) 
were consistent within 50% with the observed values (see 
Table|2]for results). In the nucleus, photoionization by stars is 
ruled out by the strong [Oi] /16300 emission ([Oi]/Ha > 0.1) 
which indicates ionization by a non-thermal source. A simulta- 
neous fit of all the line fluxes, in particular of the low [N n] A 
6583 ([Nnj/Ha' < 0.3), is obtained with the following parame- 
ters: Z/Zq = 0.5, n H = 100 cm 4 , log U = -3.4. 

In the regions located in the d,e knots, the best fit was ob- 
tained using as input ionizing source a Kurucz model atmo- 
sphere with T = 45 000 K, logn H = 2. log U = -3, Z/Zq = 
0.3. 

6.3. Photoionization + shock models 

The presence of the tidal tails strongly suggest that Mkn 298 
could be a merging system leading to collisional processes. A 
more realistic model of the emission line ratios should then 
probably take into account the coupled effect of a photoion- 
izing radiation flux from an external source and of shocks. 

j 1 

Therefore, we have adopted the code Suma ( Viegas & Contini 
1994; Contini & Viegas 2001a b, and references therein) to 
model the spectra. Suma allows to calculate the spectra emitted 
from a gas which is heated and ionized by the coupled effect of 
a photoionizing radiation flux from an external source (active 
nucleus or stars) and of shocks. The input ionizing source is a 
power-law in the AGN case. Radiation from a stellar cluster is 
used for starbursts. The input parameters are those which re- 
fer to the shock, the shock velocity, V s , the preshock density, 
no, and the preshock magnetic field, Bo, and those referring to 

3 e.g. H: 1.00, He: 0.1, C/H: 3.55 xl0~ 4 , N/H: 9.33 xl0~ 5 , S/H: 
1.62 xlO 5 , O/H: 7.41 xlO 4 



the radiation flux, the radiation intensity Fh(in units of number 
of photons cirT 2 s _1 eV~' at 1 Ryd), and the spectral index a, 
in the power-law case. In case of starbursts, the parameters re- 
ferring to radiation are the age of stellar population, t in Myr, 
and the ionization parameter, U. Moreover, the input parame- 
ters include the geometrical thickness of the clouds, D, and the 
relative abundances of He, C, N, O, Ne, Mg, Si, S, CI, Ar, and 
Fe to H. A pre-shock magnetic field Bo=10~ 4 gauss was used 
for all models. 

The shocked clouds are moving outward from the nucleus 
and radiation reaches the inner edge of the clouds opposite to 
the shock front. 

A first set of models was computed by adding a shock com- 
ponent to the active nucleus. The best fit of the observed line 
ratios was obtained with the following parameters (Table |2}: 
V s = 110 km/s, n = 100 cm" 3 , Log F H = 7.15, D = 0.36 pc. 
The abundances are nearly solar, with the exception of sulphur: 
N/H=9 x 10 4 , 0/H=6.6 x 10~ 4 ; S/H= 1.1 x 10 4 , 

In the case of a starburst + shock, the best fit was ob- 
tained with: N/H=3.5 x 10 4 , 0/H=6.3 x 10~ 4 , S/H= 5 x 10~ 6 
(Z/Zq ~ 0.6); log U = -2.52, t = 0.0 - 2.5 Myr, D = 0.5 pc, 
V s = 80 km/s, no = 60 cm 4 . Such values of the shock velocity 
and preshock dens ity are typical of t he bulk of ionized gas in 
starburst galaxies dViegas et alJl999l) . 

While in both cases the models provide a good fit of the ob- 
served line ratios, the flux from the active nucleus is extremely 
low, a factor ~ 1000 lower than what is derived e.g. in Sey fert 
2s, L og Fft > 11 or in Liners, Log Fu > 8 dContini et alJ 
2003). Moreover, the [Om]/l4363/H/3 ratio is overpredicted by 
the AGN+shock models by a factor ~ 10. 

We conclude that while the presence of an active nucleus 
cannot be excluded, it is not required in order to explain the 
observed line ratios, assuming that a shock component is also 
present. Ionization by a starburst, rather than by an active nu- 
cleus, appears to be also in better agreement with (i) the lack 
of well defined spatial trend of the line ratios with the distance 
from the putative active nucleus and (ii) the morphology of the 
nucleus derived from the images in Sec. [3] Figure|5]displays the 
distribution of the ionization parameter U and of the preshock 
density, no, which result from the best fit of the line ratios dis- 
played in Figure |8] we obtain -2.00 < log U < -2.52, 70 < 
no < 100, with higher preshock densities in Nl. This confirms 
that the starburst+shock processes are spread over the whole 
nuclear region of Mkn 298. 

The spectra of the d,e knots are all typical of starburst re- 
gions with V s =90 kms 4 , no = 50 cm 4 , log U 2.3 and an 

age of the stellar cluster t < 2.5 Myr. The best fit was obtained 
with N/H=1.5 x 10 4 , 0/H=3.3 x 10~ 4 ; S/H= 3.0 x 10~ 6 , which 
correspond to Z/Zq ~ 0.40. 

7. Conclusions 

In this paper we have presented new imaging and integral-field 
spectroscopic data analysis of Mkn 298, aimed to investigate 
its controversial nature. By means of ground-based broad-band 
R and HST band V images we have performed a detailed anal- 
ysis of the morphology of this galaxy, obtaining that Mkn 298 
is an early-type spiral, likely a SO galaxy (B/T ~ 0.78), with 
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Fig. 8. Line ratios, dust extinction and H a velocities from the rebinned MPFS spectra (R2) in the nucleus are here displayed as 
circles whose size is proportional to the measured values; the values in brackets are the minima and maxima. Overlaid is the H a 
image. 



a de Vaucouleurs bulge and without any evident nuclear point- 
like source. In addition, we have pointed out the presence of 
weak stellar shells around the galaxy. These shells, together 
with the bright and extended eastern tail, and the nuclear dust 
lane, contribute to strongly support the idea that Mkn 298 is the 
result of a major merger event. Ionized gas in Mkn 298 is dis- 
tributed in two or maybe three main knots, visible in the Ha- 
continuum subtracted image, and whose spectra show mostly 
the typical emission lines produced by thermal ionization. This 
gas is extended and oriented as the dust lane, suggesting that 
these two features are likely to be connected. The kinematic in- 
vestigation of Ha and [O ra] velocity fields show clearly that 
the gas is rotating with a maximum deprojected velocity of 



+250 km s with respect to the systemic velocity of the galaxy. 
Theref ore, the integral-fiel d data confirm the previous hypoth- 
esis bv lMoles et alJ J2004I) that a disk dominated galaxy is the 
secondary component of the merger event. 

The analysis of the emission line ratios shows that both 
AGN photoionization models and starburst+ shock models with 
sub-solar metallicity (Z/Zq < 0.6) are able to fit the ob- 
served values. In particular the high [Oi]/Ha, which could in- 
dicate the presence of an AGN, is also well reproduced by 
starburst+ shock models, where photoionization is given by a 
young (f < 2.5 Myr) stellar population and shocks are pro- 
duced in low velocity-density clouds (V s = 80 km/s, «o = 60 
cm" 3 ). AGN+shock models are ruled by the [Om]/14363/Hj6 
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ratio which is too high if compared with observed values and 
by the flux from the AGN which is unrealistically low. 

Considering that emission line ratios are constant over all 
the regions where they are detected, with no well defined spa- 
tial trend, we conclude that the starburst+shock models provide 
a more likely possibility than an AGN: merging would induce 
both collisional processes and star formation in the gas, thus 
producing the observed line ratios. No evidence was fo und for 
the very low metallicity claimed bv lMoles et alJ J2004h . Z/Zq 
< 0.1. We do not c onfirm the high SFR (~ 90 M Q /yr) found by 
iMoles et alJ J2004 : we find instead that SFR/H a < 10 M /yr. 
A comparable SFR was found taking the IRAS infrared lumi- 
nosity at the position of Mkn 298, even if this value is very 
uncertain. 

As it concerns the knots in the tidal tail (d,e), the presence 
of shocks is not so compelling as both pure photoionization 
from stars and photoionization+shock models provide a good 
fit of the observed ratios. A sub-solar metallicity (Z/Zq < 0.4) 
is found in these regions as well. 
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